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A scalar theory with a preferred reference frame is summarized. To test that theory in celestial 
mechanics, an "asymptotic" post-Newtonian (PN) scheme has been developed. This associates 
a conceptual family of self-gravitating systems with the given system, in order to have a true 
small parameter available. The resulting equations for a weakly-self-gravitating system of 
extended bodies include internal-structure effects. The internal-structure influence subsists at 
the point-particle limit — a violation of the weak equivalence principle. If one could develop an 
"asymptotic" approximation scheme in general relativity also, this could plausibly be found 
there also, in a gauge where the PN space metric would not be "conformally Euclidean" . 



1 Introduction 

The weak equivalence principle (WEP), which states that the acceleration due to gravity is 
the same for any kind of matter, applies undoubtedly to Newton's gravity. This is because 
Newton's third law implies that the gravitational self-force over a body sums to zero — as it 
may be also checked on the continuous form of Newton's attraction law, based on the Poisson 
equation. It is taken for granted that the WEP applies also to general relativity (GR), because 
the interpretation of gravitation as a purely geometric phenomenon seems to make this obvious. 
For instance, it is well-known that, by applying the local "conservation equation" for the energy- 
momentum tensor T: 

t;., v = (i) 

to dust matter, one finds that each volume element of the dust continuum has a geodesic motion. 
(More generally, Eq. (^Q) seems to imply that matter flows universally in the curved space-time, 
independently of which kind of matter is present.) However, the space-time metric, that deter- 
mines the geodesic lines, is certainly influenced by the presence of the very body whose motion is 
sought for, however small be its mass. And since this influence is nonlinear, it is not so obvious 
that it must cancel in the limit where the size of the body evanesces. Hence the motion of even 
a small dust object might after all depend, perhaps, on the structure of this object, e.g. it might 
depend on its density field and its internal motion. 

The work^ summarized here came as a natural follower of a study^ on the equations of 
motion of the mass ce nters i n a self-gravitating system of extended bodies in a scalar alternative 
theory of gravitation: BBS There, ft was found that the internal structure of the bodies does 



play a role on the motion of their mass-centers in this theory. Hence, it appeared necessary to 
investigate the point-particle limit in this theory. Moreover, with the "asymptotic" scheme of 
PN approximation (PNA)^^ which has been used in the former work, ^ the internal-structure 
influence appears rather unavoidable, quite independently of the theory considered. ^Therefore, 
in the present paper, we shall also briefly discuss what might happen in GR. 



2 Basic Principles of the "Scalar Ether-Theory" 

This is a bimetric theory, thus it endows space-time with two Lorentzian metrics: a flat metric 
7° (the "background", or "prior-geometrical" metric), and a curved metric 7 (the "physical" 
metric, which is more directly related with physical space and time measurements). It admits 
a preferred reference frame E ( "ether" ) . Most equations are written in E and are merely space- 
covariant. Space-time is assumed to be the product RxM, where M is the "space" manifold, 
i.e. the set of the positions x in the preferred reference frame E. Metric 7 is related to 7 
through the scalar field / as follows. We postulate a gravitational dilation of time standards, 
making measured time intervals smaller: 



dt^/dT = V?M = /3(x, T) < 1, (2) 

where t x is the "local time" in E, measured by a clock at a fixed point x in the frame E, and where 
T is the "absolute time" . The latter is the inertial time associated with 7 in the frame E, it is 
also the time which would be measured at a point x bound to the frame E and far enough from 
any massive body. And we postulate a gravitational contraction of material objects (including 
length standards, thus making measured distances larger than distances evaluated with the 
Euclidean metric), only in the direction of the gravity acceleration: 

J dl° ifdx = (dx l ) orthogonal tog, i.e. g°(g, dx) = 0; . . 

= I dl /y/f(x.,T) ifdx parallel tog. ^' 

Here g° is the Euclidean space metric associated with 7 in the frame E and dl° = (gr°(dx, dx)) 1 / 2 
is the corresponding line element, similarly dl is the line element defined with g, where g is the 
Riemannian space metric associated with 7, also in the frame E; and the gravity acceleration 
is: 
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S = -^of, V / = grad g o/, (grad g o/) 1 = ^fj. (4) 

(Matrix (g 0l i) is the inverse of the component matrix (goij) of the Euclidean space metric g°.) 
In coordinates bound to the frame E (and such that x° = cT), the components of metric 7 are: 
7oo = /, Hj = ~9ij, 7o« = 0, the gij 's being deduced from ©EED 

Dynamics of a test particle is governed by a relativistic extension of Newton's second law: 

DP 

F + m{v)g = —, (5) 

where Fo is the non-gravitational (e.g. electromagnetic) force, v = dx/dt x the velocity and 
v = gf(vjv) 1 / 2 its modulus; m(v) = m(0)j v is the relativistic mass (7^ is the Lorentz factor); 
P = m(v)v is the momentum; and Dw/Dt* is the appropriate time-d eriva tive of a (spatial) 
vector w in the space M endowed with the time-dependent metric g. In the static case 
(/,o = 0)> that dynamics implies Einstein's motion along geodesies of the curved space-time 
metric 7.^ For a dust continuum, we may apply © pointwise and this implies^ the following 
equation: 

t;, v = b, (6) 



(compare with Eq. Q), where 6 M is defined by 

6o(T) = i WM T^, ^(T^-i^oT *. (7) 

(Indices are raised and lowered with metric 7, unless mentioned otherwise. Semicolon means 
covariant differentiation using the Christoffel connection associated with metric 7.) Universality 
of gravity means that the same equation must remain true for any material medium. Thus, the 
dynamics of a test particle, as well as the dynamical field equation for a continuous medium, 
fully obey the "weak equivalence principle" , just as they do in GR. 

In the "non-cosmological case", the equation for the scalar field / is 

a f l (f.« \ 8lrG ,0 n 

with A = div g ograd g o the usual Laplace operator defined with the Euclidean metric g°, with G 
Newton's gravitation constant, and with a = (T 00 )e the mass-energy density in the ether frame, 
where T is the "mass tensor", i.e. the energy-momentum tensor in mass units. a 

3 Asymptotic PN approximation for relativistic theories of gravitation 

In order to obtain asymptotic expansions of the fields as function of a small field-strength param- 
eter A, one should be able to associate a (conceptual) one-parameter family (S\) of gravitating 
systems with the physically given system S. This means that one should be able to define a fam- 
ily of boundary- value problems. In relativistic theories of gravity, the natural boundary- value 
problem is the initial-value problem. This applies^ to the investigated scalar theory. Thus one 
should define a family of initial-value problems. In GR, an approximation s che me based on a 
such family of initial conditions has been introduced by Futamase & Schut2P^ although their 
initial condition for the space metric was very special. They derived (not detailed) expansions 
of the local equations but no "global" equations i. e. for the mass centers of a n N-body system. 
In the present scalar theory, detailed local equations,^ and global ones, were derived in a 
general case. 

When trying to define a family suitable for the Newtonian limit (and indeed including a 
precise definition of the latter), one is led to observe th at th ere is an exact similarity transfor- 
mation in Newtonian gravity for a perfect-fluid system: I fl l -W I If pW (pressure), = i^ 1 ) (pW ) 
(density), u^ 1 ) (velocity), and U$ (Newtonian potential), obey the Euler-Newton equations, 
then, for any A > 0, the fields 

p (A) (x,T) = \ 2 pW(x,V\T), /9 (A) (x,T) = ApW(x, V\ T), (9) 



(8) 



U$\x,T) = \U$\x,V\T), uW(x,T) = VAuW(x,\/Ar), (10) 

also obey these equations — provided the state equation for system S\ is _F( A )(p( A )) = XF^ (\~ 2 pW). 
To define the Newtonian limit in the scalar theory, we just apply the similarity transformation 

a In the most general case, applicable to a heterogeneous universe on a cosmological time-scale, the gravitational 
field is not directly the field / but instead the field of the "ether pressure" p e , and a different field equation applies. 
But this reduces to (JHJ if the time variation of the "reference pressure" p^°(T) = Sup xgM p e (x, T) is neglected.^ 
The analysis of homogeneous cosmological models^ shows that the latter variation takes place over long time 
scales indeed, of the order of 10 s years. 



to the initial data for the gravitating system of interest (with suitable modifications, in 
particular V = (c*/2)(l-/) is substituted for U n )Mt^ defines a system Sa as the solution of 
this initial-value problem Pa, and one expects that, as A — > 0, the solution fields admit expan- 
sions in powers of A, whose dominant terms have the same orders in A as in Eqs. (|9|)-(|10[). Using 
these expansions for the system of interest, S, is then justified insofar as the value Ao of A for S 
is small enough. It is then easy to check that, by adopting \M]\ = A[M] and [T]\ = [T]/vA as 
the new units for the system Sa (where [M] and [T] are the starting units of mass and time), all 
fields become order A , and the small parameter A is proportional to 1/c 2 (indeed A = (co/c) 2 , 
where cq is the velocity of light in the starting units). Hence, the derivation of PN expansions 
is easy: all fields depend on the small parameter 1/c 2 , and are thus expanded in powers of 1/c 2 , 
each with n + 1 terms (k = 0, n). Thus, each exact equation splits into n + 1 exact equations 
(it is just coefficient identification for a polynomial in A). The scalar theory admits consistent 
first-order expansions in A (or 1/c 2 ). The first term (k = 0) is Newtonian gravity, hence there 
is a correct Newtonian limit. The n = 1 approximation is the first PNA, tha t is linear in the 
PN fields (which correspond to k = 1). We call this approximation scheme^^the "asymptotic 
PNA" , because it is a rigorous application of the usual method of asymptotic expansions for a 
system of partial differential equations. 

In the standard PNA, developed for GR by Fock'^ and Chandrasekhar,^! 1/c 2 is formally 
considered as small parameter, and the matter fields p, p and u are not expanded. This implies, 
in particular, that the equations of the standard PNA are not exact ones resulting from a 
"splitting" of the starting equations into separate equations for the successive orders — as this 
is the case for the "asymptotic" PNA. In contrast with the latter, the standard PNA does not 
pertain to th e us ual method of asymptotic expansions, and the equations are different in the 
two schemes. ^El However, it does not appear easy to build a general asymptotic PNA in GR, 
for in GR the initial conditions must verify the nonlinear "constraint equations". It might be 
difficult (if it is at all possible) to build a family of initial conditions that would both satisfy the 
constraint equations and be suitable for the Newtonian limit in, say, a perfect-fluid system. 

4 PN equations of motion of the mass centers in the investigated theory 

The mass centers are defined as local barycenters of the rest-mass density, p ex act> which expands 
as 

/Wet = P{i) + 0(l/c 4 ), P(i) = p + pi/c 2 . (11) 

To get the 1PN equations of motion of the mass centers, the local 1PN equations of motion 
(which are just the first-order expansion, in A oc c~ 2 , of the continuum dynamics equation © 
as applied to a perfect fluid), are integrated inside the bodies. Integrating the zero-order and 
first-order expanded equations of the time component of the local equation of motion, in domain 
D a occupied by body (a) (a = 1, N), leads to define the 1PN mass of body (a) as 

M^=[ pmdV = M a + M^/c 2 , M a = [ pdV, M\ = [ Pl dV, (12) 

JD a JD a JD a 

(V is the Euclidean volume measure on the space M), and this integration gives just M a = 
Const., Ml = Const. Thus, the 1PN mass center a^j of body (a) is defined by 

M«a {1) = / p {1) xdV = M a a + M^/c 2 , 

with 

M a a = / px<N, M^ai = / p^aV. (14) 
JD a JD a 



(13) 



Integrating the space components of the local PN equations of motion gives the sought equations 
of motion of the PN mass center am- The 0-order equation is the Newtonian equation: 

M a a l = I pllfdV (15) 

(U is the zero-order coefficient in the expansion of V = (c 2 /2)(l - /), it obeys the Poisson 
equation with the zero-order density p; thus it is the Newtonian potential; the superscript ^ 
designates the external part of U .) The PN correction (order 1 equation) has the form 

M 1 a d\= [ f}dV. (16) 



The "PN force" density f} depends in particular on the zero-order fields: 0-order pressure p, 
density p, velocity u, Newtonian potential U. As a consequence, the internal structure of the 
bodies influences the motion already from the first PNA.^MJJ This follows naturally from using 
the "asymptotic" method of PNA and, in the author's opinion, should hold true for GR if a 
(general) asymptotic PNA could be used in GR (see the end of Sect. |3J). 



5 Point-Particle Limit. Reason for Conflict with the Weak Equivalence Principle 

We consider^the situation in which the size of one of the N bodies, say (1), is a small parameter 
£. (This body could be an asteroid or a spacecraft.) Thus we consider a family (S£) of 1PN 
systems. (Since the equations of the asymptotic 1PN approximation make a closed exact system, 
we may content ourselves with those equations and hence forget the weak-field parameter A — 
provided the system considered is indeed suitable for the PNA.) This family is defined by an 
initial data which is independent of £ apart from the size of the small body. The work consists 
then essentially in expanding the integrals on the r.h.s. of Eqs. (|15[) and (|16[) as £ — ► 0, for the 
small body, i.e. a = 1. The principal part of the Newtonian acceleration is quickly got from 
(fTTjh and (JTSJ): it is 

a = V £/ (1) (a) + O(£ 2 ). (17) 

This means that the Newtonian point-particle limit is unproblematic, as expected. To expand 
the PN correction to the acceleration, we use the simplifying assumption that the zero-order 
(Newtonian) motion of the small body is a rigid motion. 

Then it remains merely a (rather lengthy) calculus. ^ We find that, at the point-particle 
limit: £ — > 0, there remains a structure- dependent part, Ag = ord(£°), in the PN acceleration 
Am of the small body. Moreover, in the case that there is just one massive body, assumed 
static and spherical, the limit of Am differs from the PN acceleration of a test particle just 
by As- This is a patent violation of the weak equivalence principle! In the general case, the 
structure-dependent acceleration As depends on a symmetric tensor S, with 

S kj = [ pu^dN/Ah (18) 

JDi 

(where u± is the Newtonian self-potential of body (1)), and it may be written in the form 

A S = vlsmaxW, (19) 

with |w| strongly orientation-dependent, and 



A s max = 2 x 10- 7 d m/s 2 



(20) 



on Earth (d is an average density of the small body, expressed in g/cm 3 ). This magnitude seems 
dangerous, but it may not kill the theory.^ 



The general reason that makes this violation a priori possible is simply the fact that the small 
body is gravitationally active and the theory is nonlinear! Since the PN part of the coordinate 
acceleration field, iij, depends nonlinearly on (all) matter fields, we cannot even define which 
part of ui(x) comes from the body itself where x belongs, and which part comes from the other 
bodies. This nonlinear dependence is then carried over to the PN correction to the acceleration 
of the mass center, through Eq. (|14|) t. Thus, there is no a priori reason why the acceleration of 
the mass center should depend only on the other bodies in a relativistic theory of gravitation, 
even if the body is very small. Now the specific reason that makes this violation actually occur 
is the presence in the PN spatial metric of derivatives U{ of the Newtonian potential. Indeed, 
the spatial Christoffel symbols depend therefore on the second derivatives Ui j, whose self-part 
is independent on £, hence survives in the point-particle limit (see Eq. (|18j0 . But the XI i 's 
are also there in Schwarzschild's metric (this is the "anisotropy" of the latter — b y th e way, it 
is the exact solution of this scalar theory for the static spherically symm etric c ase" 1 ^). Hence, 
although the standard PN spatial metric of GR is conformally Euclidean ^ and thus no Ui 
's are there, we believe that the same violation could occur in GR, if one would (could) use 
an asymptotic PNA, and this in a gauge where the PN metric would be "Schwarzschild-like" 
instead. ^ 
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